www.gusucode.com > 用于近红外光谱中 进行校正集和验证集近红外光谱的选择源码程序MATLAB > 用于近红外光谱中 进行校正集和验证集近红外光谱的选择源码程序MATLAB/code/ks.m

    function Rank=ks(X)
%+++ Employ the K-S algorithm for selecting the representative samples;
%+++ Basic idea of KS: Find from the candidates the sample whose minimal
%    distance from the samples in representative sets is the maximal.
%+++ Hongdong Li, May 10,2008.
%+++ Revised in Dec. 5, 2009.

tic;
[Mx,Nx]=size(X);
Rank=zeros(1,Mx);           %+++ Record the indices of the representative samples.
out=1:Mx;
D=distli(X);
X=X'
[D,N]=size(X);
X2=sum(X.^2,1);
D=repmat(X2,N,1)+repmat(X2',1,N）-2*X'*X;
[i j]=find(D==max(max(D)));  
Rank(1)=i(1);               %+++ Initializes as the two samples of the furthest distance.
Rank(2)=j(1);
out([i(1) j(1)])=[];        %+++ The remaining samples.
%+++ Iteration of  K-S algorithm %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
iter=3;
while iter<=Mx
   in=Rank(find(Rank>0));
   Dsub=D(in,out);   
   [minD,indexmin]=min(Dsub);
   [maxD,indexmax]=max(minD);
   Vadd=out(indexmax);
   Rank(iter)=Vadd;
   out(find(out==Vadd))=[];
   iter=iter+1;
end
toc;
%+++ Iteration ended %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%%%%%%% END OF SUB
function D=distli(X)
X=X';
[D,N] = size(X);
X2 = sum(X.^2,1);
D = repmat(X2,N,1)+repmat(X2',1,N)-2*X'*X;